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ABSTRACT 

We present the results of deep near-infrared spectroscopy of seven submillimetre-selected 
galaxies from the SCUBA 8-mJy and CUDSS surveys. These galaxies were selected because 
they are too faint to be accessible to optical spectrographs on large telescopes. We obtain a 
spectroscopic redshift for one object, and likely redshifts for two more, based on a combina- 
tion of marginal emission line detections and the shape of the continuum. All three redshifts 
broadly agree with estimates from their radio/submm spectral energy distributions. From the 
emission line strengths of these objects, we infer star formation rates of 10-25 M Q yr _1 , while 
the lack of detections in the other objects imply even lower rates. By comparing our results 
with those of other authors, we conclude it is likely that the vast majority (more than 90 per 
cent) of the star formation in these objects is completely extinguished at rest-frame optical 
wavelengths, and the emission lines originate in a relatively unobscured region. Finally, we 
look at future prospects for making spectroscopic redshift determinations of submm galaxies. 

Key words: galaxies: distances and redshifts - galaxies: evolution - galaxies: formation - 
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1 INTRODUCTION 

The resolution of the majority of the extragalactic submillimetre 
(submm) background into individual sources (Blain et al. 1999) has 
made it possible, in principle, to study the redshift distribution of 
these sources and hence investigate the history of obscured star for- 
mation in the Universe. Such investigations are important to com- 
pare the star formation which occurs in obscured sources with that 
which takes place apparently unobscured, e.g., in the Lyman break 
galaxies identified by Steidel, Pettini & Hamilton (1995). Although 
the similar strengths of the extragalactic optical and submm back- 
grounds suggest that obscured and unobscured star formation occur 
in roughly equal quantities (Hauser et al. 1998), the more detailed 
relationship between the populations of galaxies within which these 
processes occur is unclear. Are they coeval populations, or are each 
of them the dominant mode of star formation at a particular cosmic 
epoch? 

Spectroscopy of submm- selected sources has proven difficult 
and time-consuming, and has led to the emergence of an industry to 
estimate redshifts using template SEDs which cover the wavelength 
range from the far-infrared to the radio. Since the initial work by 
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Carilli & Yun (1999, 2000), the technique has been adapted by 
Dunne, Clements & Eales (2000), Rengarajan & Takeuchi (2001), 
Yun & Carilli (2002), and Hughes et al. (2002). Often these esti- 
mates disagree strongly, however (see, e.g., Section 5 of this pa- 
per), due to differences in the assumed spectral energy distribution 
templates, and this only increases the need for spectroscopic mea- 
surements. 

There are two main issues which make it difficult to obtain 
spectroscopy of these objects. The first of these arises because the 
spatial resolution of the Submillimetre Common User Bolometer 
Array (SCUBA) is rather poor (15" FWHM), leading to a number 
of optical/infrared counterparts within the source error circle. The 
case of HDF 850.1, the brightest submm source in the Hubble Deep 
Field, serves as a good example (Hughes et al. 1998; Downes et al. 
1999; Dunlop et al. 2004). To obtain an accurate position requires 
sensitive observations with an interferometer at millimetre or cen- 
timetre wavelengths in order to achieve the arcsecond resolution 
necessary to make an unambiguous identification. 

The second difficulty arises after an accurate position has been 
determined. In many cases, the correct identification is extremely 
faint and may require deep imaging observations to be identified. 
Ivison et al. (2002; hereafter 102) report the optical/IR continuum 
characteristics of a relatively bright submm sample, having refined 
most of their positions via deep radio imaging. A large fraction was 
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shown to be sufficiently bright for spectroscopy on large telescopes, 
as later demonstrated using the LRIS-B UV/optical spectrograph 
on the Keck telescope (Chapman et al. 2003). The tendency for the 
optical/IR counterparts to be diffuse and/or multiple ("composite 
red/blue galaxies"; 102) means that the fraction of light transmit- 
ted by a spectroscopic slit is low compared to, say, a Lyman-break 
galaxy of similar magnitude. Continuum magnitudes are not a re- 
liable indicator of likely success, however, if targets display strong 
emission lines, as many have been shown to do, including several 
with / > 26 (we use Vega magnitudes throughout). Nevertheless, 
the ability to detect continuum emission with a reasonable signal- 
to-noise ratio is essential if emission lines are weak or lie outside 
the available wavelength coverage, and conclusions about the un- 
derlying stellar population can only be drawn if there is sufficient 
signal-to-noise ratio in the continuum to detect breaks and/or stellar 
absorption feautures. Given the optical faintness and frequent red 
colours of these sources, sensitive near-infrared spectroscopy may 
prove to be a fruitful avenue for sources where optical spectroscopy 
has proven unsuccessful. 

In this paper, we present the results from deep near-IR spec- 
troscopy of a sample of submm galaxies, made with the extremely 
sensitive OHS instrument on Subaru Telescope. The format is as 
follows. In Section 2, we describe the selection of our sample and 
our observations, and in Section 3, we detail the careful data reduc- 
tion procedures required to extract maximum sensitivity from these 
faint objects. Section 4 presents our results, while in Section 5 we 
describe the spectra of individual objects in detail. We discuss our 
results in Section 6, and speculate on future prospects for obtaining 
complete spectroscopic redshifts for submillimetre samples in Sec- 
tion 7. Finally, we summarize our results and conclusions in Sec- 
tion 8. Throughout this paper, we adopt Ho = 70 krns" 1 Mpc -1 , 
On = 0.3, and ft A = 0.7. 



2 SAMPLE SELECTION AND OBSERVATIONS 

Our sample of targets was taken from the SCUBA 8-mJy Survey 
of the Lockman Hole East and ELAIS N2 fields (Scott et al. 2002; 
Fox et al. 2002) and the Canada-UK Deep Submillimetre Survey 
(CUDSS) 14-hour field (Eales et al. 2000; Webb et al. 2003). We 
selected only targets which had been securely identified through 
follow-up radio observations which were able to provide sufficient 
astrometric accuracy to isolate a single object 1 . Furthermore, we 
selected targets for which the identification had very red optical/IR 
colours (see 102), so as to exclude galaxies for which optical spec- 
troscopy might be expected to be more successful. The sample is 
therefore by no means complete, and was constructed to test the 
limits of the most sensitive near-infrared spectrograph currently 
available on any telescope and assess the likely success of future 
IR spectroscopy. 

Galaxies were observed with the OH-airglow suppressor 
(OHS; Iwamuro et al. 2001) and the Cooled Infrared Spectrograph 
and Camera for OHS (CISCO; Motohara et al. 2002) on Subaru 
Telescope during the period UT 2002 May 20-24. OHS is a spec- 
troscopic filter which eliminates 224 strong OH airglow lines in the 
J and H bands from a 1.3 x 28 arcsec 2 region of sky and outputs 
an undispersed beam which is then fed to CISCO to enable im- 
ages or spectra to be taken. The broad-band sky background in this 

1 ELAIS N2 850.12 was not formally detected in the radio, though 102 note 
that the combination of its red colour and the presence of several 3cr radio 
peaks makes the identification reasonably secure. 



Table 1. Observing log. The dates refer to the spectroscopic observations. 
Imaging observations were made immediately prior to the spectroscopy, 
except for LE 850.1, which was imaged on UT 2002 May 2 1 (see Section 
5.1). 



PA Exposure (s) 



Name 


UT Date 


O 


Imaging 


Spectroscopy 


LE 850.1 


2002 May 23 


167 


12 x 50 


8 X 1000 


LE 850.3 


2002 May 20 


43 


6 x 20 


8 x 1000 


CUDSS 14.1 


2002 May 20 


90 


6 x 30 


8 x 1000 


CUDSS 14.3 


2002 May 24 


118 


6 x 50 


8 x 900 


CUDSS 14.9 


2002 May 21 


63 


6 x 20 


8 x 900 


ELAIS N2 850.2 


2002 May 23 


30 


6 x 50 


8 x 1000 


ELAIS N2 850.12 


2002 May 24 


82 


6 x 50 


8 x 900 



suppressed region (the 'dark lane') is reduced by a factor of ~ 25, 
enabling the acquisition of extremely faint targets. Conditions were 
clear for all the observations presented here, and the optical seeing 
was measured to be 0.5-0.7 arcsec from the auto-guider camera. 

Each target was observed in an identical manner. First, a 
nearby offset star was acquired and centred in the OH-suppressed 
'dark lane', using CISCO in imaging mode with an i/-band filter. 
The telescope was then moved according to the offset between this 
star and the target, as measured from existing near-IR images. The 
location of the target on the array was then determined by nodding 
the telescope along the slit and taking images at each position, the 
number and length of which depended on the expected magnitude 
of the target. The images were combined to remove the sky back- 
ground and allow the target to be seen. If the target was not within 
the dark lane, an appropriate telescope offset was applied and the 
process repeated until an accurate centroid could be determined. A 
small offset was then applied using the autoguider camera to move 
the target into the centre of the dark lane. The spectrograph slit was 
then closed to 1 arcsec and aligned with the dark lane, and the JH 
grism inserted into the light path. Eight exposures were then taken 
in two quads using the standard 'ABBA' technique (e.g. Eales & 
Rawlings 1993), with exposure times of 900 or 1000 seconds per 
position. The offset between positions was usually 10 arcsec, al- 
though we chose to align the slit to include any nearby objects and 
a smaller offset was used if the the standard value would cause con- 
fusion between the beams. In some cases, a small offset along the 
slit was applied between the two quads. When these exposures had 
finished, the telescope was slewed to a nearby atmospheric ratioing 
standard which was acquired without opening the slit by making 
small telescope movements perpendicular to the slit position angle 
until the count rate from the star was sufficiently high. The stars 
used were HIP 52256, HIP 71078, and HIP 80419 (Perryman et al. 
1997) for the Lockman Hole, CUDSS 14-h, and ELAIS N2 fields, 
respectively. Eight exposures (of 5-10 seconds each) were taken in 
a single quad with two frames per position. 

Observations were also made of argon and neon calibration 
lamps and the planetary nebula NGC 7027, to provide wavelength 
calibration and to measure the spectral resolution of the instrument. 
These observations were made with a 0.5-arcsec slit to more allow a 
more accurate determination of the locations of the emission lines. 
No shifts were found in the location of any spectral features, or in 
the locations of the unsuppressed night-sky lines, during the course 
of the run. An observing log is presented in TableQ 
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3 DATA REDUCTION 

The low level of signal coming from the extremely faint targets 
(we stress that our targets are producing detector signals of about 1 
electron per pixel per minute, which is less than the detector dark 
current) makes it imperative to remove all detector artefacts. To do 
this, we first subtracted one image from each pair from the other 
to remove the contribution from the dark current. Next, any differ- 
ences in bias level between the four quadrants of the detector were 
removed by determining the signal in unexposed regions of each 
quadrant and subtracting this constant value from all the pixels in 
that quadrant. Variance images for each combined pair were con- 
structed based on the raw pixel values in the original images and 
the detector read noise. Residual sky lines were removed by fit- 
ting a low-order function to the columns of the images, excluding 
any regions contaminated by objects. The two pairs in each quad 
were then combined, by a straight average if the pixel values in 
each pair were consistent at the 5u level, or by the value closest to 
zero if they were not. Next, the image was shifted by an amount 
equal to the offset between the positive and negative beams and 
this shifted image was subtracted from the original, with further re- 
jection of deviant pixels. The variance images were combined in 
a statistically-appropriate manner. Then the two quads were com- 
bined and spectra of all objects extracted in 1-arcsec (9-pixel) aper- 
tures. The pixel-by-pixel uncertainty on this spectrum was deter- 
mined from the variance image, although this is likely to underes- 
timate the true noise level in regions of strong sky emission where 
the sky level may not be accurately determined. 

Spectra of the atmospheric ratioing standards were obtained 
by simply averaging all the spectral images taken in the 'A' posi- 
tion and subtracting the average of the spectral images taken in the 
'B' position. This image was shifted to align the two beams and 
subtracted from the original version of itself. A spectrum was ex- 
tracted in a 1-arcsec aperture. 

A quadratic wavelength solution was determined from the arc 
and planetary nebula spectra and applied to all the spectra. The 
r.m.s. deviation of the individual emission lines from the fit was 
~ 9 A (approximately 1 pixel). Model spectra of the standard stars 
were constructed assuming they could be described by blackbod- 
ies with effective temperature and _ff-band magnitudes determined 
from the K-band magnitudes and spectral types listed in the HIP- 
PARCOS catalogue (Tokunaga 2000). Spectra of the targets in flux 
units were then determined by dividing each target spectrum by the 
appropriate atmospheric standard spectrum and multiplying by the 
model blackbody. Since the standard stars were chosen to be close 
in the sky to the targets, their airmasses at the time of observa- 
tion are similar to those of the target at the end of its observation, 
rather than an average target airmass. This may result in imperfect 
removal of atmospheric absorption features. 

The absolute flux calibration of the targets was determined 
from the acuqisition imaging. Images at each of the two posi- 
tions were combined and then subtracted from each other. This 
image was then shifted to align the two beams and subtracted 
from the original version of itself. Photometry was performed in 
a 1-arcsec square aperture and therefore the flux scale of the final 
spectra represents the amount of light in the extraction aperture. 
The uncertainty in the photometry was estimated by fitting sepa- 
rate Gaussians to sigma-clipped, logarithmically-scaled histogram 
of the pixel values both inside and outside the dark lane. The instru- 
mental zeropoint of the photometry was tied to observations of the 
UKIRT Faint Standards FS 27 and FS 146 (Hawarden et al. 2001), 
taken at the end of the last night. Any systematic uncertainties in 



the absolute calibration introduced by only taking such images at 
one point during the run are likely to be dwarfed by the 10-20 per 
cent random uncertainties in the target photometry. In cases where 
more than one object was in the spectrograph slit, the flux scale was 
tied to the brightest object. 

Inevitably, we expect some spurious features from 'warm' pix- 
els which are not rejected by our rejection algorithm, since this 
algorithm needs to be conservative to avoid throwing out a large 
amount of good data (in addition, the rejection limits need to be 
loose enough to not flag a change in the signal which is the re- 
sult of changes in the seeing and/or atmospheric transparency). The 
reduced two-dimensional spectra are therefore used to determine 
the veracity of any putative emission lines, by determining whether 
they arise from only one of the two beams. 



4 RESULTS 

The FWHM of the arc lamp and planetary nebula lines was mea- 
sured to be 35-40 A, indicating that the object spectra will have a 
spectral resolution of R ~ 200. This is consistent with directly- 
made measurements of the resolution obtained with a 1-arcsec slit 
(e.g. Willott et al. 2003), and means that any emission lines will 
be unresolved. The low resolution of the spectrum allows us to re- 
ject as spurious some features which might at first glance appear 
real (especially in the smoothed spectra presented in Fig.0. This 
is important because the large number of pixels means that we have 
a ~ 50 per cent chance of finding a 3<j 'emission line' simply by 
chance. This probability is reduced by a factor of about 2 by exclud- 
ing those 'lines' where the flux is dominated by a small number of 
errant pixels. On the other hand, it is almost inevitable that each 
spectrum will contain at least one emission feature which is for- 
mally significant at the 2.5a level, and we must treat such marginal 
'detections' with extreme caution. 

We can therefore only consider redshifts to be plausible if they 
are determined from a single, highly significant line, or from two 
or more less significant detections whose relative wavelengths cor- 
respond to a plausible pair of features. To detect such features, we 
use a 51-pixel running median to estimate the continuum level, and 
then search for groups of 9 pixels (i.e. ~ 75 A, or approximately 
the resolution FWHM) whose total deviation from this continuum 
has a signal-to-noise ratio larger than a certain threshold, which was 
set at 2.2 after some experimentation. Once features had been de- 
tected in this manner, the 2D spectrum was examined to ensure that 
they were not the result of poor sky subtraction or low-level cosmic 
rays which had not been removed by the rejection algorithm. Not 
all detections are listed therefore in Table|2| where we summarize 
the results of our analysis, although we only exclude those features 
which are clearly spurious. The detections are described in detail 
on an object-by-object basis in the next section. It is important to 
note that the uncertainties quoted for the line fluxes are derived 
solely from the detector read noise and photon shot noise, and do 
not take into consideration the effects of improper sky subtractor or 
the placement of the continuum. Incorrectly determining the con- 
tinuum by as little as 0.1 x 10 -17 W m -2 fim -1 (see Fig.Q) will 
cause the measured line flux to be in error by ~ 10~ 20 Wm -2 , 
which is typically 20-40 per cent of the flux. We therefore warn 
against any interpretation of line ratios, which are often substan- 
tially different from those expected. Even though in some cases the 
median appears by eye to be wrong (and will always be slightly bi- 
ased by the presence of the detected line), any 'by eye' estimate is 
likely to be just as wrong, so we choose to fix the continuum at the 
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Figure 1. Spectra of the targets, extracted in 1 arcsecond wide apertures. 
The solid line is the target spectrum, smoothed with a 5-pixel boxcar filter, 
while the dotted line indicates the lcr uncertainty in the smoothed spectrum. 
The greyscale at the top indicates the OH-suppressed sky brightness - re- 
gions of high sky brightness are subject to additional noise from imperfect 
sky subtraction, which is not accounted for in the plotted noise spectrum. 
The brightest sky features (black in the greyscale) are multiplet O2 tran- 
sitions, which are not rejected by the mask. Emission lines are labelled in 
those objects for which we have derived a redshift. 
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Figure 1. continued 



level derived by the running median, which at least has the advan- 
tage of being non-subjective. The line centres are determined from 
a simple centroiding of those pixels above the continuum level, and 
typically have an accuracy of ~ 10-20 A which, when combined 
with the residuals from the wavelength solution fit, imply uncer- 
tainties of ±0.005 in each individual redshift measurement. 

Given the marginal nature of our line detections, we investi- 
gate the probability of obtaining false redshifts from our data. We 
first produce 100 artificial spectra from each of the ten noise spec- 
tra (seven objects, plus three companions), by multiplying the noise 
spectrum by a spectrum of random numbers drawn from a Normal 
distribution. We apply our line detection algorithm to each of these 
spectra, again rejecting those detections which are not consistent 
with unresolved lines. Using the 2.2a detection threshold, we find 
that most spectra produce one or two spurious features, but there 
is a pronounced tail, and 3 per cent of our spectra produce five fea- 
tures. While this may seem like cause for alarm, the fact that we 
require more than half of the detected features to correspond to 
plausible lines (and be absorption or emission features as appropri- 
ate) makes it extremely unlikely that we will derive a redshift from 
five randomly-located lines. 

We investigate two distinct ways in which we can derive an 
erroneous redshift. First, the redshift can be produced from entirely 
spurious lines. This is not possible unless there are at least two 
lines, and in this case the redshifts derived from their identifications 
must agree to within 0.015 (for the case of three spurious lines, we 
require that the redshifts from any pair agree, while for the cases 
of four or five spurious lines, we require that consistent redshifts 
are derived from a triplet). These probabilities can be calculated di- 
rectly from our 1000 spectra, and amount to 10 per cent. Second, 
an erroneous redshift can be produced from a single real line, plus 
one or more spurious lines. To calculate this probability requires a 
priori knowledge of the probability of there being a genuine fea- 
ture at a given wavelength. If we assume that the genuine emission 
line has an equal probability of appearing at any point in the spec- 
trum, then the probability of an erroneous redshift is 20 per cent, 
multiplied by the a priori probability that there is a genuine line in 
the J or H bands. The most likely cause of an erroneous redshift 
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Table 2. Results of our line detection algorithm. 
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Figure 2. Spectra of 'companion' objects, presented in an identical manner 
to the targets of Fig.fTI 



arises from cases where three lines are detected, and a pair of them 
correspond to plausible features. 







Flux 






Object 




(10~ 20 Wm~ 2 ) 


Identification 


2 


LE 850.1 


1.1544 


4.3 ±1.8 








1.2440 


2.0 ±0.9 








1.5998 


1.8 ±0.8 






LE 850.3 


1.1627 


3.0 ±1.4 


[O II] A3727 


2.120 




1.5993 


1.6 ±0.7 


spurious 




CUDSS 14.1 


1.2994 


1.9 ±0.8 








1 6675 


-1.6 ±0.6 


See text 






1 7227 


2.3 ±1.0 


for discussion 




1.7543 


2.1±0.7 






cunss 14 3 




Nothing detected 




CUDSS 14.9 


1.5216 


-1.6 ±0.6 


CaK? 


2.869 




1.5374 


-2.0 ±0.8 


Ca H? 


2.874 


N2 850.2 


1.2889 


4.0 ±1.8 


[O II] A3727 


2.458 




1.5637 


-2.4 ±0.8 


spurious 






1.5751 


2.6 ±1.2 


spurious 






1.6777 


2.1 ±0.9 


H/3 


2.451 




1.7091 


3.4 ±1.2 


[O III] A4959 


2.446 




1.7304 


3.9 ±1.3 


[O III] A5007 


2.456 


N2 850.12 


1.5722 


2.0 ±0.8 


spurious 






1.6678 


2.2 ±0.8 




2.431 




1.7125 


2.5 ±0.8 


[O III] A5007 


2.420 


Companion objects 


CUDSS 14.3 




Nothing detected 




N2 850.2 


1.1700 


3.8 ±1.6 







N2 850.12 



Nothing detected 



5 NOTES ON INDIVIDUAL OBJECTS 

In this section, we describe our redshift determinations based on 
the data in Table|2| We list all plausible emission lines, even if they 
do not provide a definitive redshift measurement, so that our data 
may be of use if an independent redshift determination supports the 
reality of these lines. We deliberately refrain from making subjec- 
tive comments on the relative believability of any lines in a given 
spectrum. 



5.1 LE 850.1 

Due to the extremely faint and diffuse nature of this source, it took 
a substantial amount of imaging on the night of UT 2002 May 2 1 
to make a firm imaging detection. Although by the time we had 
done so, the hour angle of the source was no longer favourable, we 
confirmed the offsets from a nearby star as being 11.6 arcsec west 
and 5.2 arcsec south. We returned to LE 850. 1 two nights later, and 
used our previous imaging to hasten the acquisition. Since the dis- 
tance between the target and offset star is larger than the OHS field 
of view, we first located the bright star near the end of the field of 
view, and made a small (2 arcsec) movement of the telescope by al- 
tering the location of the centroiding box on the autoguider camera. 



The success of this technique depends on the accuracy to which the 
plate scales of CISCO and the AG camera are known, and we esti- 
mate the error to be no more than 0.2 arcsec. Since we were unable 
to obtain reliable photometry from our acquisition image, we ap- 
plied the same flux scaling factor for LE 850. 1 as was applied for 
LE 850.3; from the flux in the 'clean' part of our ff-band spectrum 
(1.60-1.70 (im), we estimate H = 23.1 ± 0.4 in our extraction 
aperture. 

Three emission lines are suggested by our line-finding algo- 
rithm, although all are extremely marginal. No pair of these lines 
corresponds to plausible features, and so we must conclude that at 
least two, and possible all three, of the lines are spurious. However, 
given that it is three times more likely for there to be only two spu- 
rious lines, rather than three, we consider it probably that one of 
the lines is real. The redshift estimate of z — 2.61q \ (model 'le2' 
of Aretxaga et al. 2003; hereafter A03) suggests that either of the 
bluest two lines could be [O II] (at z = 2.097 or z = 2.338), but 
we are unable to assign a redshift from our data. 

5.2 LE 850.3 

This source is the brightest of our targets in the ff-band (H = 
20.4±0.2 in our extraction aperture) and has relatively good signal- 
to-noise in the clean parts of the spectrum. Our detection algorithm 



6 C. Simpson et al. 




4000 4500 5000 5500 



Figure 3. Spectrum of LE 850.3 (as Fig. Q, overlaid (dashed line) with 
the spectrum of a 250Myr-old starburst with an e-folding time of 100 Myr, 
reddened by Ay = 2.2 mag. 

produces two emission lines, both at marginal (< 2.5<r) signifi- 
cance. The lines do not correspond to any plausible pair of tran- 
sitions, but we note the presence of a clear continuum break at 
~ 1.2 fim. We confirm the reality of this break by subtracting the 
spectrum of blank sky at a distance of 1.5 arcsec along the slit (ex- 
tracted in an identical 1 -arcsec aperture). If we identify this break as 
the Balmer jump, then the bluer of the two emission lines is [O II] 
at z = 2.120 ± 0.004, consistent with the radio-submm estimate 
of 102). 

Fig.[3]shows the near-IR spectrum of LE 850.3 (shifted to the 
rest frame), overlaid with a synthetic starburst model (from the PE- 
GASE code of Fioc & Rocca-Volmerange 1997). This model has an 
decaying star-formation rate with an e-folding time of 100 Myr and 
is seen at an age of 250 Myr. It has been reddened by Av = 2.2, 
according to Pei's (1992) parametrization of the Milky Way extinc- 
tion law. Although this is not a fit, it provides a reasonable descrip- 
tion of the observed continuum. The presence of a strong Balmer 
jump requires that the rest-frame 4000-A light be dominated by 
~ 10 8 yr-old stars, although this may be because younger stars are 
more heavily obscured (rather than because the star-formation rate 
is decaying with time, as in our simple model). The reasonable fit 
that this model provides to the spectrum gives us confidence in as- 
signing a redshift of 2 = 2. 12 to this source, even though it is rather 
inconsistent with the estimate of z = 3.0t_a'g °f A03. 

5.3 CUDSS 14.1 

This source (also known as CUDSS 14A) was the subject of in- 
tense study by Gear et al. (2000). High-resolution Hubble Space 
Telescope and ground-based (Canada-France-Hawaii Telescope) 
adaptive optics images reveal it to be a red, compact galaxy. Deep 
near-IR JHK spectroscopy by Gear et al. failed to detect any emis- 
sion lines to a level of 1-2 x 10 -19 W m -2 . Our new spectroscopy 
pushes this limit almost an order of magnitude fainter, but fails to 
find any significant (> 3<r) lines. We detect four features (three 
emission lines, and one absorption line) of marginal significance, 
however. 

The measured wavelengths of the bluest and reddest emis- 
sion lines are consistent with being H/3 and Ha, respectively, at 
z = 1.673. In this interpretation, the third emission line and the 
absorption line are artefacts. Alternatively, these two lines could 
be [O 11] at z = 2.486 and [O in] A5007 at z = 2.504; al- 
though this redshift discrepancy is large, it is not implausible given 
the low signal-to-noise ratios of the detections. With these iden- 
tifications, the third emission line could be A4959 at z — 2.474, 
which is barely consistent with the redshift from the other two 
lines, while the absorption feature is, again, spurious. The lower 
redshift is consistent with both the photometric redshift estimate 



of z = 1.25 ± 0.3 and the radio-submm redshift estimate of 
z = 1.9 ± 0.48 (Clements et al. 2004; hereafter C04). We note, 
however, that A03 suggest a much higher redshift of z = 3.81q's> 
which is just consistent (at the 90% level) with both the above sets 
of line identifications. If this redshift is correct, then one of the H- 
band emission lines could be [O II] at z ~ 3.7, while the remain- 
ing three features would be artefacts (possibly the 1.2994-/im line 
could be Mg II A2798 at z = 3.644 if the 1.7227-/X line is [O 11] 
at z = 3.622, but Mg II is not expected unless there is an AGN 
present). Ultimately, we decide against assigning a spectroscopic 
redshift to this object. 

5.4 CUDSS 14.3 

Apart from a spurious peak due to poor subtraction of the unsup- 
pressed OH lines in the range 1.50-1.52 /im, there are no signif- 
icant features in the spectrum of this faint (H — 21.6 ± 0.3) 
object. A brighter, slightly bluer object is located approximately 
2 arcsec west and slightly north of CUDSS 14.3 (fig. 4 of Webb 
et al. 2003), whose spectrum is similarly featureless (Fig.[2j. We 
are unable to obtain a redshift for either object. We note that both 
the photometric and radio-submm redshift estimates of C04 are 
z = 1.1 ± 0.3, at which redshift no strong emission lines would 
lie within our spectral coverage (Ha would fall between the J and 
H spectroscopic windows). A03 obtain a much higher redshift es- 
timate (z = 3.2~^I q), which would place the [O II] line some- 
where between the 1.27-^tm O2 night sky emission and the long- 
wavelength end of our spectral coverage, and it could therefore eas- 
ily fall between the J and H windows, or near the blue end of the 
//-band where the poor sky subtraction affects our sensitivity. 

5.5 CUDSS 14.9 

The spectrum of CUDSS 14.9 (H = 20.9±0.2) bears aremarkable 
similarity to that of LE 850.3, with an apparent break in the con- 
tinuum around 1.2 /im. However, the spectrum is noisier and the 
break is weaker - if we make an ad hoc definition od the break 
strength as B = f, ' 2 , Mm S\ dA/ f, 1 ' 1 f Mm S\ dA, we measure 
B = 1.08 ± 0.37 for CUDSS 14.9, compared to B = 1.47 ± 0.24 
for LE 850.3 (or 1.85 ± 0.39 after removal of the [O II] line). 

Two lines are detected by our algorithm; absorption features at 
1.5216 and 1.5374 fim with observed equivalent widths of ~ 40 A 
and ~ 50 A, respectively. The wavelengths of these lines are con- 
sistent with the Ca H+K lines at z « 2.87. We are hesitant about 
assigning this redshift, however, since there is no clear continuum 
break associated with the absorption lines (although there is insuf- 
ficient coverage to the blue of them to make an accurate measure- 
ment of the continuum level) and the redder of the two features 
is affected by a night sky emission line. Furthermore, this redshift 
is inconsistent with C04's photometric and radio-submm estimates 
of z = 1.44 ± 0.32 and z = 1.7 ± 0.43, respectively, although 
it agrees with A03's estimate of z — 2.51j'q. If our spectroscopic 
measurement is correct, however, none of the strong rest-frame op- 
tical emission lines associated with star formation (Ha, H/3, [O II]) 
would lie within a spectroscopic window. 

5.6 ELAIS N2 850.2 

The identification for this source is extremely faint (H = 22.4±0.4 
in our spectroscopic aperture) and the continuum is only seen in 
our 2D spectrum after substantial smoothing. The location of the 
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extraction apertures were determined from the acquistion image. 
Five plausible emission features (one of which is detected at > 3cr) 
and one absorption feature are detected by our algorithm, although 
the absorption feature is close to one of the emission lines and may 
result from imperfect subtraction of the 1.58-^im sky emission. The 
other four emission lines, however, are consistent with a redshift 
z — 2.453 ± 0.006. This redshift is consistent with tha lack of 
formally-detected continuum in the J-band, since the 4000-A break 
will lie between the J and H atmospheric windows, and also with 
the redshift estimates of 102 and A03. It also agrees well with the 
optical redshift of 2.454 (Chapman et al., in preparation). 

Our spectroscopic slit was aligned to pass through a galaxy 
2arcsec south-southeast of ELAIS N2 850.2 (fig. 2 of 102). This 
source shows faint continuum in the clean part of the J-band 
(1.20 < A < 1.26 (im; the continuum is more clearly seen in 
the smoothed 2D spectrum than in Fig.|2j but this does not rule 
out its being at the same redshift as the submm source, since it 
is brighter (H = 22.1 ± 0.3 in our aperture). There is a single 
marginally-detected emission line at 1.1700 /im; if real, this would 
mean that the 'companion' is not physically associated with the 
submm source, since it does not correspond to any known emission 
line at the redshift of the submm source. 

5.7 ELAIS N2 850.12 

This identification for this source is also very faint (H = 22.2 ± 
0.4) and the signal-to-noise ratio of our spectrum is low. However, 
we find three emission lines, of which we identify two as H/3 and 
[O III] A5007 at z = 2.425 ± 0.005 (the absence of detectable 
[O III] A4959 is simply due to its low flux). This redshift is also 
consistent with the lack of continuum at the red end of the J-band, 
since it places the 4000-A break at 1.37 /im (we note, however, 
that continuum is detected at the ~ 3u level in the 'clean' part 
of the J-band (1.20-1.25 urn). At this redshift, the [O II] A3727 
emission line would lie at 1.276 ± 0.001 /im, where our sensi- 
tivity is low (3<r » 7 x 10" 20 Wm" 2 ) due to the 2 night-sky 
emission; its non-detection is therefore not at odds with a ratio 
[O Il]/H/3 ~ 3 (e.g. Kennicutt 1998), especially if there is some 
reddening. We note that this redshift is unexpectedly low, given 
the high ratio of submm to radio flux, though this inconsistency 
can be understood if the handful of 3a radio peaks noted by 102 
are the relics of a diffuse radio source. Since this redshift is de- 
rived from two of three detected lines, it should be treated with 
caution, although it gains some support from the fact that there is 
positive flux at the expected position of [O II] (formally, the value 
is (3.8 ±2.2) x 10~ 20 W m~ 2 ), and from the relatively high sig- 
nificance of the line at 1.7125 /im (S/N > 3). 

Our spectroscopic slit was aligned to pass through a second 
object located approximately 7arcsec west of N2 850.12 (on the 
edge of the postage stamp shown in fig. 2 of 102). This object is 
brighter (H = 20.2 ± 0.1) and bluer and our algorithm failed to 
detect any emission lines in its spectrum (see Fig. [5). We cannot 
say whether it is physically associated with the submm galaxy. 



6 ON THE ABSENCE OF EMISSION LINES 

Given the extreme star-formation rates implied by the submm lu- 
minosities of these objects and our extensive wavelength cover- 
age, the lack of strong emission lines is perhaps a surprise. A star- 
formation rate of 1000 Mq yr _1 (typical of these objects; Ivison et 
al. 1998; Gear et al. 2000) corresponds to an Ha luminosity Lho, = 
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Figure 4. Sensitivity to emission lines associated with star formation. The 
bottom line indicates our typical 3cr sensitivity to an unresolved emission 
line as a function of wavelength. The three curves near the top of the plot 
show the wavelengths and line fluxes expected for a 1000 Mq yr — 1 star- 
burst (Kennicutt 1992, 1998). Vertical tick marks are placed at redshift in- 
tervals of 0.1, and are numbered at intervals of 0.5 Finally, the downward- 
pointing vertical arrows indicate the effect of Ay = 3 mag of internal 
reddening. 

1.3 x 10 37 W (Kennicutt 1998), and we assume Ho/H/3 = 2.78 
(Case B; e.g. Osterbrock 1989) and Ha/[0 II] = 0.9 (Kennicutt 
1992). These lines should all be well above our detection threshold 
(see Fig.|4j. In fact, in the absence of any extinction, our 2-hour 
integrations should be sensitive to star-formation rates as low as 
~ 20 Mq yr~ out to z ~ 3.5, although there is a gap in our cov- 
erage between 2.6 < z < 3.0 when [O II] lies in the gap between the 
J and H atmospheric windows. This gap could be significant since 
it lies close to the peak in the redshift distribution, z mc dian ~ 2.4, 
presented by Chapman et al. (2003), based on optical spectroscopy 
of bluer radio-detected submm galaxies. 

We have redshifts for only three of our seven targets. Taken 
at face value, the measured line fluxes of LE 850.3, N2 850.2, and 
N2 850.12 correspond to star-formation rates of ~ 10, ~ 20, and 
~ 25MQyr _1 , respectively, well below the estimates from the 
radio-submm data (102). Such large discrepancies are not unusual; 
for example, in the extremely red object HR 10 (Hu & Ridgway 
1994), the raw (i.e. uncorrected for extinction) Ha-derived star- 
formation rate is a factor of ~ 30 lower than that derived from the 
submm flux (Cimatti et al. 1998; Dey et al. 1999). The two rates can 
be reconciled by appealing to substantial quantities of dust which 
severely attenuate the observed line emisison. In all three cases, the 
emission-line-derived star-formation rate can be rectified with the 
submm-radio-derived SFR if the emission lines are extinguished 
by Ay ~ 3-4 mag (compared to the Ay m 4.5 mag needed in 
HR 10 2 ). Slit losses are unlikely to account for more than a factor 
of about 2 in our line fluxes since our slit width corresponds to a 
linear size of 8 kpc, which is large for a starburst region (Smail et 
al. 2003). Furthermore, the radio emission is not observed to be 
extended on scales much larger than this. 

2 Dey et al. (1999) claim that a simple foreground dust screen is inappro- 
riate because it would imply an intrinsic ratio of [O II]/Hq 1.3, which 
they claim is "at the extreme limit of observed values for local star-forming 
galaxies (e.g., Kennicutt 1992)." However, Fig. 12 of Kennicutt (1992) in- 
dicates that this ratio is quite normal, since the data plotted are uncorrected 
for reddening (a correction factor of ~ 2) or [N II] contamination (~ 1.5). 
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Arguing against a foreground screen hypothesis is the detec- 
tion of Lya emission in submillimetre galaxies (Chapman et al. 
2003). The measured fluxes correspond to star formation rates of 
~ 10 M Q yr~\ broadly in line with our measurements, while a 
foreground screen with Av = 3 would produce at least four or- 
ders of magnitude of extinction at Lya and completely suppress 
the line. The consistency between the UV and optical star forma- 
tion diagnostic lines therefore suggests that we are seeing a small 
fraction (~ 1 per cent) of the star formation in a relatively unob- 
scured way, while the vast majority is completely extinguished at 
optical/UV wavelengths. Some obscuration is still required since a 
10 M yr _1 SFR at z = 2.4 produces an observed V w 22 (Lei- 
therer, Robert, & Heckman 1995), which is brighter than our targets 
(102; C04). This picture is quite reasonable, since massive stars are 
formed in the local Universe inside the cores of dense molecular 
clouds, and the ~ 15 Myr it takes for these stars to disperse their 
parent clouds is longer than the lifetimes of stars with M > 15- 
20 Mo (see discussion in Jimenez et al. 2000). Jimenez et al. sug- 
gest that the observed flux of far-ultraviolet photons represents only 
about one-sixth of the total star formation, and the observed line 
emission should also underestimate the true star formation rate by 
a similar factor if the massive stars are heavily obscured. Larger dis- 
crepancies are possible if the starburst is younger than 15 Myr since 
even the longer-lived stars are embedded in their parent clouds. 

An alternative hypothesis which warrants investigation is that 
the star formation rates derived from the submm observations are 
substantially overestimated. The most likely reason for this would 
be erroneously high assumed dust temperatures. For a dust emis- 
sivity index f3, the bolometric luminosity per unit dust mass is pro- 
portional to T 4+/3 , while the dust mass for a given monochromatic 
submm luminosity is inversely proportional to the dust tempera- 
ture. Therefore, for a source with a known redshift and measured 
850-/im flux, the inferred bolometric dust luminosity (and hence 
star formation rate) is proportional to T 3+l3 (e.g. Blain et al. 2002). 
This suggests that very large errors could be introduced by an in- 
correct estimate of the dust temperature. However, the expression 
is only accurate at low redshift, where hv c <C kT {v c being the 
frequency at which the observed 850-pim radiation is emitted), and 
the possible error is more modest at z ~ 2, where a factor of 2 
uncertainty in dust temperature corresponds to a factor of 10 uncer- 
tainty in luminosity (e.g. Hughes, Dunlop & Rawlings 1997; Eales 
et al. 2000). 

102 adopt T dust = 45 K, which is near the hot end of the dis- 
tribution found for ultraluminous infrared galaxies (Dunne et al. 
2000; Farrah et al. 2002). In principle, therefore, the star forma- 
tion rates of ~ 1000 Mq yr -1 could be an order of magnitude too 
high if the temperature of the dust is ~ 20 K. The ratio of 850-/im 
to 450-/im flux is a relatively sensitive function of dust tempera- 
ture for z > 1, as the shorter wavelength point approaches the peak 
of the emission. Eales et al. (2000) claim an average upper limit of 
Ssso /S450 ^0.5 for the CUDSS sources, which appears to rule out 
T ~ 40 K dust if the median source redshift is 2. However, Fox et 
al. (2002) only constrain this ratio to be > 0.3 for sources from the 
8-mJy survey, which is low enough to allow hotter dust. In prin- 
ciple, the radio flux should provide an alternative estimate of the 
SFR if the redshift is known, since the tight radio-FIR correlation 
provides a way to derive the total FIR luminosity without having to 
assume the dust temperature. Systematically overestimated SFRs 
should result in high values of S&50 M m/<Si.4 ghz (i.e. steeper spec- 
tral indices) whereas there appears to be a strong tendency for the 
observations to produce shallower spectra (e.g. C04). The cause of 
this is unclear, and it seems that this issue will only be unequivo- 



cally resolved with shorter wavelength detections which can con- 
strain the shape of the SED near the peak of the thermal emission. 
For now, we believe we are justified in claiming that it is highly 
improbable that the submm-derived SFRs are reconcilable with our 
spectroscopically-derived estimates unless most of the star forma- 
tion is heavily obscured. We are unable to determine whether this 
obscured fraction is 90 per cent or 99 per cent of the total star for- 
mation. 



7 COMPARISON WITH OTHER SPECTROSCOPIC 
STUDIES AND FUTURE PROSPECTS 

Near-infrared spectroscopy of submm galaxies has been rather lim- 
ited to date. Ivison et al. (2000) presented A"-band spectra of two 
sources behind the lensing cluster Abell 1835, SMM J1401 1+0252 
and SMM J14010+0253, which both showed Ha emission but 
lacked the strong [N II] emission characteristic of an AGN. Us- 
ing near-infrared spectrographs on larger telescopes, Frayer et al. 
(2003) presented a spectrum of SMM J04431+0210, while Smail 
et al. (2003a,b) observed SMM J17142+5016 (in the protocluster 
around LBDS 53W002) and ELAIS N2 850.4, respectively. All 
three of these galaxies showed strong emission lines indicative of 
the presence of an AGN, which Smail et al. suggest is a common 
feature of these objects. All the emission lines in these objects are 
brighter than our upper limits, by an order of magnitude in the cases 
of SMM J17142+5016 and N2 850.4. The spectra of our targets 
are dramatically different from those of the (few) objects in the lit- 
erature. However, this fact may be readily explicable by the fact 
that our targets were selected to be optically faint, and are there- 
fore perhaps biased towards more heavily extinguished objects. On 
the other hand, previously observed objects already had known red- 
shifts from optical spectroscopy, and must therefore have had bright 
Lya emission, and/or sufficiently bright optical continua to allow 
the detection of absorption lines. Either case implies that the ob- 
jects with near-IR spectra in the literature are lightly extinguished. 

In support of this idea, Eisenhauer et al. (2003) recently ob- 
tained a high-quality near-infrared spectrum of the z = 2.565 
submm source SMM J1401 1+0252. They found that the rest-frame 
optical continuum of this object is dominated by stars with ages 
of a few hundred Myr, similar to our estimate for LE 850.3 (the 
only object in which we detect an age-sensitive stellar feature) al- 
though, unlike LE 850.3, the continuum of SMM J1401 1+0252 ap- 
pears to be unreddened, which presumably explains its unusually 
bright rest-frame UV spectrum (Ivison et al. 2001). Such an age is 
perhaps larger than expected, since a constant star formation rate of 
~ 1000 Mq yr _1 can form a massive galaxy in only ~ 100 Myr. 

The most serious difficulty which arises when attempting to 
understand the nature of submm galaxies as a class comes from a 
lack of clarity in the criteria with which submm galaxies are se- 
lected for spectroscopy by different groups. This difficulty is al- 
most certainly compounded by the likelihood that spectra are only 
published where firm redshifts are obtained (i.e., objects with lumi- 
nous emission lines, or those with bright enough continua to allow 
the detection of absorption features). This leads to the presence of 
strong and uncertain selection effects in the literature samples, and 
makes generalized statements about the properties of the popula- 
tion as a whole unreliable. 

Nevertheless, almost 100 spectroscopic redshifts for submm 
galaxies have been obtained, and there is little reason to believe that 
the redshift distribution obtained from these objects is significantly 
different from that of the class as a whole (Chapman et al. 2003; but 
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Figure 5. Simulated 7-hour FMOS spectrum of ELAIS N2 850.2. The spike 
at 1.62 fim is caused by the masking of multiple airglow lines which effec- 
tively remove all the flux from that pixel, and is not significant (as indicated 
by the dotted line which shows the lcr noise). See text for more details. 

see also Dannerbauer et al. 2004). It is therefore clearly time for the 
study of these important objects to progress from simple redshift 
determination and basic classification (AGN, starbust, composite) 
to more detailed analyses. The difficulties posed by the presence of 
large quantities of dust and the existence of multiple components 
point to the need for spectroscopy over a long wavelength baseline, 
covering rest-frame wavelengths from Lya to Ha. 

The most interesting new instrument for this purpose is 
FMOS, the Fibre Multi-Object Spectrograph (Maihara et al. 2000; 
Kimura et al. 2003) which is scheduled to be commissioned on 
Subaru Telescope in late 2005. This instrument will provide OH- 
suppressed spectroscopy with increased sensitivity, spectral reso- 
lution, and wavelength coverage (0.9-1.8 /im) compared to OHS, 
using 400 fibres across a 30-arcminute field of view (this corre- 
sponds to the source density at a flux level of Ssbo ^ 3mJy; e.g., 
Borys et al. 2003). The ability to take many spectra simultaneously 
enables much longer exposures than we were able to undertake here 
and, coupled with the increased sensitivity of the new instrument, 
it should be relatively straightforward to probe a further order of 
magnitude below our current detection limits. Our emission line 
sensitivities will therefore correspond to star formation rates of 
~ lMoyr - , and it would be remarkable if a significant num- 
ber of targets failed to display emission lines. Furthermore, the in- 
creased sensitivity to the continuum emission will allow more ac- 
curate modelling of the galaxies' stellar populations. As a demon- 
stration of the likely data quality we can expect from FMOS, we 
simulated the spectrum of N2 850.2 using the FMOS Spectral Sim- 
ulator developed by Naoyuki Tamura 3 (Fig[5J- Our input spectrum 
assumes fluxes for the H/3, [O II], and [O III] A5007 lines given 
in Table|2| superimposed on a continuum of a 3-Gyr-old elliptical, 
reddened by E(B — V) = 0.5. All emission lines are assumed 
to have FWHM = 500 kms -1 , and an exposure time of 7hours 
was adopted, representing the maximum likely on-source integra- 
tion time possible in a single night's observing. 

A major caveat to this lies in the infrared "spectroscopic 

3 The FMOS Spectrum Simulator can be found at 
http://elvira.phyaig.dur.ac.uk/naoyuki.tamura/simulator.html 



desert" at 2.6 < z < 3.0. Here, [O II], the Balmer jump, and the 
4000-A break are redshifted between the J and H windows; H/3 
and [O III] are redshifted between H and K; and Ha is redshifted 
beyond the K band. Although this is a relatively small redshift 
range, it is close to the peak of the redshift distribution and may 
include more than 10 per cent of sources from a flux-limited sam- 
ple (Chapman et al. 2003). Without the reddening-insensitive age 
measurements provided by the Balmer jump and/or 4000-A break, 
and with the low spectral resolution of FMOS, detailed stellar pop- 
ulation analysis of these targets is unlikely to prove fruitful. 

For now, further infrared spectroscopy is still worthwhile as 
a precursor to the above observations. Many submm sources dis- 
play optical emission lines and, regardless of whether these are 
powered by a starburst or an AGN, there will be associated neb- 
ular lines redshifted into the near-infrared. The relative strengths of 
these lines (in objects whose redshifts are suitable for near-infrared 
spectroscopy) will test our idea that they come from a relatively un- 
obscured region comprising only a small fraction of the total star- 
burst. 



8 SUMMARY 

We have presented the results of deep near-IR spectroscopy of 
seven galaxies selected from the SCUBA 8-mJy and CUDSS 
submm surveys. We have obtained spectroscopic redshifts for three 
of these objects, two of which (LE 850.3 and ELAIS N2 850.2) 
agree well with the prior estimates from radio-submm data, while 
the other (ELAIS N2 850.12) is unexpectedly low. These redshifts 
are based on marginal ( < 3er) detections of emission lines, al- 
though in all cases they come from more than one spectral feature, 
and in one instance (N2 850.2) there are four emission lines, and the 
derived redshift agrees with the (as yet unpublished) Lya and CO 
redshifts. In all three cases, the star-formation rates derived from 
the emission-line luminosities are much lower than those derived 
from the radio-submm data. By comparing the bulk properties of 
our sample with those observed at optical wavelengths, we suggest 
that we are seeing a small fraction of the total starburst through low 
extinction. Further infrared spectroscopy of sources with secure op- 
tical redshifts are required to test this hypothesis. 
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